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We study the dynamics of a particle in continuous time and space, the displacement of which 

is governed by an internal degree of freedom (spin). In one definite limit, the so-called quantum 

random walk is recovered but, although quite simple, the model possesses a rich variety of dynamics 

and goes far beyond this problem. Generally speaking, our framework can describe the motion of an 

electron in a magnetic sea near the Fermi level when linearisation of the dispersion law is possible, 

coupled to a transverse magnetic field. Quite unexpected behaviours are obtained. In particular, 

we find that when the initial wave packet is fully localized in space, the Jz angular momentum 

component is frozen; this is an interesting example of an observable which, although it is not a 

^^ ', constant of motion, has a constant expectation value. For a non-completely localized wave packet, 

r— ^ . the effect still occurs although less pronounced, and the spin keeps for ever memory of its initial 

f^ ' state. Generally speaking, as time goes on, the spatial density profile looks rather complex, as a 

fvj ^ consequence of the competition between drift and precession, and displays various shapes according 

to the ratio between the Larmor period and the characteristic time of flight. The density profile 

gradually changes from a multimodal quickly moving distribution when the scatttering rate is small, 

to a unimodal standing but flattening distribution in the opposite case. 
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I. INTRODUCTION 



> 

in 

^^ ' The term quantum random walk has been coined some times ago by Aharonov and Davidovichi to qualify the 

^^O , motion of a quantum particle moving either left or right according to the value of the Jz component of its spin. As 

^-P such, this is completely different of the quantum Brownian motion problem, in which one fully describes in a quantum 

framework the motion of a single (light) particle moving in a bath of heavy particles - a quantum version of the 

classical Brownian motion. This latter problem has been the subject of many papers in the past, revivified in the 80's 

i-G [ by the pioneering work of Caldeira and Leggett (for a review, see Leggett eta&), showing that for the so-called ohmic 

coupling with the bath, dynamical symmetry breaking occurs above a given threshold. On a more elementary level, 

simple models show how the classical brownian motion can be generalized to the quantum world, although long time 

scales (and algebraic behaviours) appear in correlation functions at small enough temperature, as shown by Aslangul 

^ etal^^. We shall adopt here the recent meaning as proposed by Aharonov and Davidovich, shortened as QRW in the 

qh following. 

JV ' This latter problem is presently the subject of an intense activity, following the basic papers by Ambainis et 

.,—f [ a&, Nayak and Vishwanath^ and Diir et a* (for a recent comprehensive review, see the preprint by Kempe^). The 
^^ 1 basic motivation lies in the fact that QRW yields a dynamics at the opposite of classical random walk: for a spin 
?H I 2 particle, and a symmetric initial spin state with definite phases (see below for details), the probability density in 
space Pqrw(2;, t) displays two well-defined peaks, moving away linearly in time from the starting place. This is to 
be contrasted with the t^^^ spreading of the classical packet, the particle staying (in the average) at its starting point 
when no drift is present. It can be said that QRW, except for oscillations in spatial density, displays the same basic 
features as the classical non-difffusive motion with a drift (which can be obtained from the biased diffusion equation 
^ = —V[vP — DV P] - in the (singular) limit D ^ 0, as explained for instance in the book by Gardiner—) -, when 
the velocity itself is a random variable assuming two values ±iio with definite probabilities. This rather intuitive 
picture has been firmly grounded by the recent work of Blanchard and Hongler— . 

Due to the quantum nature of the walk, all steps are strongly correlated (as contrasted to classical motion) , as 
it is also the case for repeated mesurements yielding the Zeno paradox theoretically found by Misra and SudarshaniS 
and observed by Itano et a^}^\ this means that specific ever- lasting correlations are always relevant. One consequence 
is that the linear dimension of the visited space increases linearly in time, i. e. much more rapidly for QRW than for 
classical motion. It is hoped, on a somewhat speculative level, that Monte-Carlo algorithms could be improved by 
drawing benefit of this fact (see e. g. Kempe^^) . 
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We adopt here a more general point of view. Indeed, we define a simple model containing a single parameter, 
noted a in the following, which is essentially the product of the spin-flip rate by the time of flight in space. As 
shown below, the limit a — > can be viewed as the continuous space-time version of the conventional QRW. Out 
this limiting case, the model can describe the dynamics of an electron near the Fermi level (once the dispersion law 
has been linearized, e{k) — > zLh{k — kp)vF as is done e. g. in the Luttinger model^'^), moving in a magnetic sea 
created by impurities, or coupled to nuclear spins (Vagneni^). Elastic collisions with the magnetic background can 
flip the spin of the moving electron and, simultaneously, change its velocity from +vf to —v-p. Alternatively, such 
an reversal can be induced by a transverse magnetic field forcing the spin to precess harmonically at the Larmor 
frequency w. As a whole, spin and orbital degrees of freedom are coupled, yielding intrication of the state as time goes 
on, and competition between precession and translation in space. This produces interesting effects; the first one (and 
probably the most unexpected) is freezing of the precession when the initial packet is quite narrow - an effect which 
could have interesting applications in spintronics, as well in two-dimensional systems (McGuire et ai^), in carbon 
nanotubes (Yang et a'M.), in quantum dots (Levitov and Rashbail^) and in semiconductors (DyakonoAiiSi) . A second 
characteristic feature is the oc t spreading of the wavepacket which arises in any case, even when the precession is 
extremely rapid: because of the latter, the particle has no time to choose between right and left moves, but dispersion 
still occurs, and very quickly as compared to classical random motion. 

Obviously, the physical relevance of our model is subjected to small deviations around fcp. On the other side, 
it is hoped that the basic (and surprising) results found below keep some relevance even with such restrictions, and 
at least can take place on space and time scales to be precised in real life. It is worthy to note that, provided such 
a physical picture is meaningful, a becomes an easily tunable parameter by varying cither the magnetic field, or the 
filling-ratio of the band, or the density of the magnetic diffusing centers. 

This paper is organized as follows. We first define the model (section ^J, we explicitely show the relation 
between our model and QRW and we give a qualitative discussion of limiting cases. We then briefiy derive the 
equations giving the time-evolution of the density for an arbitrary spin J f section |lll|l . on which the spin- locking due 
to space confinement in the a — > limit can be directly shown. This fact is confirmed by a detailed study of the 
average value J^ fsection llV|) . Then, we focus on the J = 5 case (sectionEJi and calculate the density profile P{x, t), 
which displays many various shapes, some unexpected when the two time-scales related to the Larmor precession and 
the displacement in space arc of the same order of magnitude. Eventually, conclusions are drawn and hints for future 
work are given. In the Appendix, details are given on the asymptotic analytical work required to get insight on the 
features of the probability density. 

II. MODEL AND QUALITATIVE DISCUSSSION 

For a spin J particle, our model Hamiltonian reads: 

H = Lojy + rr^v.~pj,_ . (1) 

In this expression, lu is the spin-flip rate due to scattering on impurities, or the Larmor frequency due to a transverse 
magnetic field along the y direction, wis the (scalar) quantity defining the velocity-scale, p is the momentum (p = —iKS/ 
in the ^-representation) . Despite some similitude at first glance, the Hamiltonian (Q has nothing to do with the Dirac 
Hamiltonian, in which the momentum is coupled to the Dirac a matrices; indeed, the angular momentum (spin) is 
there not given by a and, in addition, the problem has here the dimensionality 2 J -|- 1 in spin space, instead of four 
in the Dirac theory. Also note that, in its one-dimensional form, with a space dependent uj and for J — 1^, this 
Hamiltonian has been used in various contexts such as inhomogeneous supraconductors (de Gennea^Si) and solitons in 
polyacetylene (Takayama et aSS.). 

In the following, we restrict to one-dimensional space; calling Ox the line on which the particle moves, the 
Hamiltonian Q simplifies to: 

d 
H = ujJy + h~ vpJz , p ~ — iTiT^ • (2) 

ox 

Note that the label x of direction of motion and the three directions defining the components Jx , Jy and Jj generally 
have no relations between them. Also note that the "kinetic" term is invariant under time-reversal. 

A. Continuous limit of the Quantum Brownian Walk 

In order to perform the continuous limit of QRW, let us recall the basic formalism. In the original model^, the 
spin-i particle hops on a lattice (lattice spacing a, n € Z) from one site to the two first-neighbors, either right or 



left according to the value + or — of the J^ component. In obvious notations, the operator S generating this spatial 
motion is 

S= \+){+\®Y.\^+l){n\ + \-){-\®Y.\^-l){n\ . (3) 

n n 

Once a jump is done, the spin is changed by the action of a Hadamard matrix T acting on the |±) spin states, which 
we choose of the following form: 



V2 



1 -1 
1 1 



(4) 



Thus, one step of the motion is generated by the product TS^ which reflects the basic sequential nature of the walk. 
As a whole, for integer times t, the state |^(i)) obeys the following equation: 

\^{t+l)) - TS\^{t)) = ifQRw|*(t)> . (5) 

Performing now a Fourier transformation in space (|fc) — N^^^^ J2n e''^""), one obtains: 

Hqkw = ^e-"'"5^«e-"'"^'^^^|fe)(fc| . (6) 

k 

Let us now formally replace t + I hy t + At, ^ by ujAt in the first exponential; taking the limit a ^ 0, At ^ 0, 
-^ ~ Const = V, one has: 

i^QRW -> 1 + —Y.i'^Jy + l^vJzMikl . (7) 

fc 

Going back to direct space, and taking now the limit liniAt^o sj [|*(3^7 ^ + ^0) ~ 1^(2;, t))], © yields: 

in^^\^{x,t)) = H\^ix,t)) (8) 

where H is the Hamiltonian given in |(2Jl. This Schrodinger equation retains the two essential features of QRW: the 
direction of the motion is determined by the value of the Jz component, and the latter is not a constant of motion 
due to the transverse magnetic field (external or intrinsic) to which the particle is coupled though the operator Jy. 

Yet, a basic difference exists between -ffqRw and H, due to the fact that in the discrete version, the particle 
jumps, and then the spin is changed by T: as already mentionned, the rules of the game are essentially sequential, 
allowing to state that the spin changes slowly as compared to the time of flight. By contrast, with H given in |(2Jl, 
both motion and spin-flip occur simultaneously. This means that, within the framework defined by H, one expects to 
recover the quantum random motion only when the Larmor frequency lo is quite small as compared to the time scale 
of the displacement. As it stands, the Hamiltonian H defines a model for which ordinary QRW is just one limit. 

One additional ingredient is required, namely the initial state, which will be chosen of the spin-space separate 
form: 

+,/ 

M=-J 

\M) is the eigenstate of Jz with the eigenvalue Mh; for physical purpose, ■0(x) is chosen as a localized even function 
with a width a, assuming real values in order to avoid any built-in inessential drift; for explicit calculations, we retain 
the gaussian normalized form: 

iPix) = (%/2^a)-i/2e-^'/(^'"') . (10) 

Obviousy, the state becomes intricate as time goes on, and one generally has at time t > 0: 

+.1 

\^ix, t)) = J2 v^A^(^' *)i^^) ■ (11) 

M=-,7 



The main goal is to find the probabiHty density in space, given by: 

+./ 

The model is now completely defined, and involves just one dimensionless parameter, called a in all the following: 

a = — . (13) 

V 

When a is large, the spin undergoes many Larmor precessions during a relatively small displacement in space; on the 
contrary, a small means that the spin precesses quite slowly when moving in space. 

B. The limiting cases 

Let us now briefly describe the limiting cases. As explained above, the a ^ limit must reproduce QRW. 
More precisely, when the Larmor period becomes much larger that any other relevant time-scale, the limit of the 
spatial density is the (2 J + l)-modal distribution: 

+J 
lim P(x, t) = y \cm\^ \^{x ~ Mvt)\^ \/t (14) 

a — > -^ — ^ 

(for J = i, the two-peak splitting of ordinary QRW is recovered). Such a density trivially gives {x){t) = 
vt^j^jhr (Jj)(0) and Ax^(i) = {vt)'^h~ AJ|(0). The ex t^ increase of its mean-square deviation merely reflects 
the fact that each peak of the density moves away linearly in time due to the persisting multimodal character of 
the density profile, which is frankly different from spreading in the usual sense. Note that, in the limit a = 0, the 
relative phases of the cm play no role, a symmetry which is broken for any finite a: in the general case, this phases 
are essential (see below) and, in particular, determine whether the density P{x, t) is symmetric in space or not. 

On the other side, when a goes to infinity, spin and space degrees of freedom become decoupled. Thus the spin 
stays immobile, and simply precesses within its initial wave packet, which remains as it stood at the beginning: 

lim P{x, t) = \ip{x)\'^ Vt . (15) 

OL — * -t-OO 

It will be seen in the following that these two trivial limits are indeed singular, especially the limit a ^ (because 
the velocity v is in factor of the highest derivative in H)\ this can be considered as a first symptom of the richness of 
the dynamics for any finite a. Anyway, the above limiting behaviours of P{x, t) are expected to hold approximately 
true in the general case for times i «C ^21 and t <^ ^, respectively and mimick the actual dynamics. 

III. FORMAL EXPRESSION OF THE DENSITY 

The Schrodinger equation is formally easily solved by going to the p-representation. It then reads: 

i;i^|$(p, t)) = [iojy + h-'vpj,) mp, t)) , (16) 

where |$(p, t)) is the p-representation of the state at time t [p is now a scalar). The time-evolution operator U{p, t) 
is such that: 

|$(p, t)) = U{p, i)|$(p, 0)) , (17) 

where \^{p, 0)) = (/)(p) |x) is the ^-representation of the initial state 0; with the gaussian (|1U|) . one has: 



By introducing the p-dependent unitary transformation R(j)) — e*^'''^ s{p)j^ ^ U{p, t) assumes the form: 

U{p,t) = i?t(p)e-'^"'"(P)*^-i?(p) , (19) 



where: 

nij,) = |u^ + (vh-^pf]i , (20) 

This allows to write down the formal expression of the Fourier transform of the density probability, V{k, t) = 
J_^ dxe^'^^P{x, t), as the following: 

dp0;</._(x|t/t(p+— ,t)[/(p-— ,t)|x) , (22) 

where (j)± = 0(p ± hk/2). For arbitrary J, such an expression seems untractable as it stands, but it allows to look 
at the limiting case a — > 0, which can be obtained by assuming a fully localized packet (a = 0+). Starting with the 
initial gaussian wavepacket IjlUII . a careful limiting procedure yields: 



;i/;;d,.-.-(* 



BmP(t, f) = ,/- / d(e-'i ixW'"' ■'■U) ■ (23) 

From this, one immediately obtains the limiting expression of the probability density in direct space: 

+J 
liia P{x,t)= Y^ \cm\'^ S{x - Mvt) . (24) 

"^ M=-J 

This says that the initial narrow packet splits off in 2 J + 1 components, each of them going away from the starting 
point with its own velocity Mv. Expression (|24|l holds true for any J and any initial spin state. 

This result is at first surprising (it would also trivially occur in the absence of magnetic background or of the 
transverse field {oj = 0) - which also gives a = 0). It means that extreme confinement {a = 0) of the spin forces the 
component J^ to have a constant expectation value, although Jz is not a constant of motion due to the fact that u> ^ 0. 
Obviously enough, one can question the validity of the above limiting procedure, because the value a = is indeed 
singular; in fact, the above result, which provides an oversimplified picture of the spin- freezing phenomenon, can be 
easily confirmed by analysing an innocent- looking observable (indeed Jz itself), directly obtained by the Heisenberg 
equations. This is done in the following section. 

IV. DYNAMICS OF THE COORDINATE AND OF THE COMPONENT J^ 

In order to get more insight in the above result, we now solve the Heisenberg equations of motion. They write: 

Xii = vh,~^Jzli , ph = , JzH^ -i^JxH ■ (25) 

JxH = ^JzH-fl~^VpHJyH , JyH ^ fl^^VpaJxH ■ (26) 

These equations can be readily integrated to yield: 

Jzait) ^ e{t)-J , (27) 

XH{t)=x + vn-'^f{t).J , (28) 

where J (resp. x) coincides with Jh(0) (resp. a;H(0) and where the components of the scalar vectors e(i) and T(t) 
are: 

^x (t) = — sin 9 sin Qt , 

ey{t) ^ sine cos 6(1 -cos nt) , (29) 

Cz {t) = cos^ + sin^ 6 cos fit , (30) 



and Tu{t) — Jg e„(i') dt' . The mean-square deviation of the coordinate is: 

Ax^it) = Ax2(0) + (^;?^-l)'^((T„T„J„J„) - {TuJu){T,J,)) 



(31) 



where the brackets denote averages over the initial state (see © and p8|) ') - remember that 17 and are functions of 
p, see (|20|l and (|21|l . and enter in a convolution with the Fourier transform of l|10() . For any separate initial state, the 
averages factorize: (T„J^) = {Ty}{Jy) and so on. 

These results allow a straightforward discussion displaying the strange features of the dynamics, especially the 
rather counterintuitive role of the initial spin state on the subsequent motion, especially on the symmetry of the 
probability density, as already discussed (see e. g. the analysis by Kempe^). 

Close inspection first shows the essential role of the phases of the coefficients cm appearing in the expansion 
0. Indeed, it is readily seen that when the initial spin state is an eigenvector of Jy, then {Jz)(t) and {x)(t) are 
constant (and thus vanish at any time), whereas Ax^ is still oc t^. Another consequence is that the density probability 
P{x, t) is symmetric in space in this case, and only in this case. Thus, the dynamics, and the parity of the spatial 
density, strongly depend on the relative phases of the coefficients cm defining the initial spin state (see Q), not only 
of the weights |cmP - a feature which clearly separates the general a-finite case from the a = 0+ limit. For any other 
preparation, P(— x, t) ^ -P(x, t) and the expectation values {Jz){t) and (x)(i) do vary in time, as examplified below 
for a definite preparation. 

Indeed, let us assume that the initial spin state is the eigenstate \M) of Jz {Jz\M) — Mfi\M))] then, the 
expectation values at time t are: 



{J^){t) = Mh{cos'^0 + sm^0cosnt) , 



(32) 



sin fit , 



(x) (t) ^ Mv t(cos^ 9 + sin^ 0—p^—) 



nt 



(33) 



_[j(j + 1) _ M^]{T^ + T^) + M\{T^) 



In the preceding equations, (•) denotes the average over the orbital variable 



(Tz)') 



Mj--- 



dp . 



(34) 



(35) 




FIG. 1: Variation of the time average r; (see (I37l l) as a function of a. 
In particular, the expectation value of Jz is given by {Jz){t) — ri{t){Jz){0) with: 



+00 g-2C 



2 , ^2 > ^' + "' COS Va' +e— dC 



■ Ujt 

a 



(36) 



(note that, with our definitions, — = — ). r]{t) (obviously bounded by ±1) is clearly an oscillating function of time. 
This expression allows to discuss the unexpected behaviour of the average value of Jz when a varies. First, let us 
look at the time-average value oi rj^rj = limt^+oo j /q ri{t')dt' . One readily finds: 



where $ is the probability integralSi. rj has the following behaviours: 
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FIG. 2: Variation in time of r;(i) for a narrow initial wave paclcet (left) and for a large one (right). For an infinitely narrow 
wave packet, ( Jz) remains constant. The unit of time is one period T of the Larmor precession. 

In addition, it is readily seen that 1 — 2?^ < ii{t) < 1. Thus, for a small, T]{t) oscillates around a value which is 
quite close to one, showing that (J^) becomes nearly constant in time. On the contrary, for a large, the oscillation 
takes place symmetrically around a quite small value. The variation of ^ as a function of a is given in Fig. ^ 

At short times, one has ri{t) ~ 1 — ^(ujt)'^. The behaviour of j]{t) at large times is easily found by using a 
saddle-point method. We find: 



rjit) ~ ?7 + 



2a 



'(-*+l) 



. U — 1 \ 

, t 3> mm ( — , OJ ) 

V 



(39) 



The enveloppe of the oscillation around the time averaged value rj thus decreases as i^^/^. The variation of ri{t) at 
any time is plotted in Fig|21for two values of a. 

These results confirm the confinement-locking of the spin: whereas Jz is not a constant of motion, narrowing 
the width of the initial wave packet yields an expectation value which is less and less varying in time. 

From (P5|) . the average position of the particle is easily calculated and can be written as: 



{x){t) = Mvt 



a 



-t-oo 



-2C 



(q,2 + ^2)3/2 



sin V Q;^ 



a 



(40) 



This represents a drift motion, with a damped- in-time (oc t~^) oscillation around the central position Mvtrj, imaging 
the forward/backward motion of the particle within its wavepacket as the spin precesses, causing inversion of the 
velocity. 

It is readily seen that (x) {t) is bounded for any time and any a; one finds: 



Mvt^- [uoty^] < {x){t) < Mvt[Ti+ {ujty^] 



(41) 
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FIG. 3: Variation in time of {x){t) for four values of a, when the initial state is an eigenvector of Jz. Note the different 
horizontal and vertical scales. 

This shows that, for such a preparation, the motion is always ballistic and simply approaches {x}{t) ~ Mrjvt at times 
which are large compared to the precession time scale. However, note that the effective velocity rjv decreases rapidly 
when a becomes large (see (j38|l and Fig. |21l. In the asymptotic regime, one precisely has: 



{x){t) ~ Mvt 



2a . I TTN 



(42) 



Note that in the limit a ^ but for an arbitrary initial spin state, one has {Jz){t) — {Jz){0) and: 

{x){t)=vth-\j^){0) , Ax^{t)=v^t^h-'^AJl{Q) , (43) 

in agreement with H24() . The ex t^ increase of the mean-square dispersion merely reflects the fact that in this limit, 
the density is just the superposition of the 2 J + 1 5{x — Mvt) functions. 

Obviously, the mean square deviation gives a first feeling about the spatial density, but the qualitative discussion 
given above convinces one that a given t^ increase of Aa;^(i) can be realized in a variety of ways, ranging from two 
moving sharp peaks to a gaussian-like flattening in situ. Clearly, a more precise analysis of the profile is required, and 
this is done in the following section for the J = \ case. 

V. TIME EVOLUTION OF THE SPATIAL DENSITY IN THE SPIN i CASE 



In the J = \ case, the evolution operator can be easily written explicitely. After some algebra, we find the 
propagator U{p, t) as the following: 



n(p)t n(p)t ^, , VL(p)t ^, , 
U{p, t) — cos 1 — 1 sm cos 0(j)) cr^ — i sm sm Oijp) Uy 

where the cr„'s are the Pauli matrices. From this, one readily obtains the amplitudes i-'±{x, t), given by: 

1 



(44) 



■)p±{x, t) 



'2Tih J^ 



dpe 



.px 



n{p)t . n{p)t\ ■ nr , ■ ^ip)t 

cos — - — =F 1 cos 0{p) sm — - — c± =F sm 0[p) sm — - — Czf 



^ip) 



(45) 



Close inspection of the integrals in (|45|l again reveals that for any initial state which is even and real, the density 



P(x, t) = J2e=± \''p£i^j ^)l is even only when |c+| 



id -^ purely imaginary, i. e. when the initial spin state 



is an eigenvector of Jy. In all other cases, P{—x, t) ^ P{x, t). 
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FIG. 4: Snapshots of the (symmetric) density profile at successive times t for various values of the parameter a; each curve 
is labelled by -^ where T is the Larmor period ^. The abscissa unit is the width a of the initial gaussian packet. Note the 
different scales from one drawing to another. 

In the following, we restrict to the symmetric case, taking c^ = -y= and c_ = ^, corresponding to |x) = l+)a- 
Then, the density can be written as follows: 



P(X, t) = -^ Y. \^r{x, t)P , 



(46) 



where the three quantities Ir are given integrals. With the gaussian wave packet H18|) . the latter explicitely write 



h{x, t) = ^ / dC ^== sin y^^Ml^- 



(47) 
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(48) 
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(49) 



These expressions allow an easy numerical calculation of the density P{x, t) in various cases; the results are shown 
in Fig. ^ and display the extreme variety of P{x, t) when the single parameter a varies. For a small (remember 



10 



this corresponds to QRW), the (two) peaks structure is clearly visible, and displays small satellites at the back of the 
packet. They are easily understood as arising from a (quantum) path in which the particle has undergone a small 
number of precessions. This is confirmed by the ct — ^ curves: it can be checked that the number of peaks is close 
to 4, where T is the Larmor period. The extreme-right peak arises when the particle did not precess at all, then 
comes a secondary peak associated with one precession, and so on. When a increases, the structure is still present 
but diminishes quickly as time goes on. Eventually, for a large enough, the profile does not display any structure and 
vaguely looks like a standing gaussian wave packet. As a whole, the profile is rather sensitive to a; also note that 
P{x, t) remains notably non-zero in the central region even at large times (see below, especially (|52|) '). 

Approximate analytical expressions of the density can also be obtained (see Appendix). For a <C 1, we find 
that: 



P(x, t) 



16(27r)3/2 



+ 00 



dX'c-i(^-^'-f)'i/f'(X') 



(50) 



where X' = aJ^\X'\ and where H\ is the conventional Hankel function. Details given in the Appendix allow to 

understand that in this case (a <C 1), the two main peaks moving with the velocity | are accompanied by small 
satellites, corresponding to those quantum paths with a few precessions. 

On the other side, for a 3> 1, one finds the plain gaussian distribution: 



P(x, t) 



2A(t)^ 



27rA(i) 



A = fj 



/ w^i 



\Aa'^uj 



1/2 



(51) 



Because the precession is rapid, the particle does not move in the mean, but the effect is not averaged to zero and 
produces a linear in time increase of the width of the distribution. 

Generally speaking, it appears that, for any a, the central region remains populated due to the fact that 
P{x = 0, t) decreases rather slowly in time. Indeed, a stationnary-phase argument shows that: 



P{x = 0, t) 



1 



27rA(i) 



V 



V t > Max 



27r a 

LO ' V 



(52) 



As for any initial state, the mean square dispersion of the coordinate always increases ex i at large times, and 
can be readily calculated for the symmetric case, using the results of section Hvl In particular, one finds: 



lim 

t — ^+oo 



Aa;(t) 



cos^ 0)1/2^ 



V(a)v 



with: 



(cos 



a 



27ra(3 + 4a2)e2" [1 - $(y2a)] 



(53) 
(54) 



V(a) 




FIG. 5: Variation as a function of a of the velocity V of the standard deviation Ax(i) (see 1551 1 ) 
V{a) is a monotonically decreasing function of a (see Fig. O, with the following behaviours: 



V{a) 



i(l-3yfa2) a<l 



\/3 



a> 1 



(55) 
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This means that, apart from small damped in time oscillations, the motion is essentially ballistic, with a velocity 
V{a) going to zero when the precession frequency increases. With the reduced units used in Fig^ the width of the 
distribution is 2TT/{aV{a)) (t/T); rough estimates show that the gross features of the density profile are in agreement 
with the latter expression for the standard deviation. 




10 20 30 40 50 60 



FIG. 6: Variation in time of the Shannon entropy S{t) (left). The right part shows the ratio S{t)/ In i, which becomes a constant 
at large times. 

As a final measure of the profile, let us analyse the Shannon entropy, defined as: 



/ + 00 
dxP{x, t)lnP{x, t) 
-OO 



(56) 



Three typical plots are given in the left part of Fig. |^ showing that S(t) - an ever-increasing function in time -, 
changes gently with a and does not display a transition reflecting the bimodal/unimodal cross-over. In addition, it is 
seen that S{t) oc In t at large times, which means that, as it is often the case, S{t) ex In Aa;(t). 



VI. CONCLUSIONS 



In this paper, we presented a simple model describing the dynamics of a particle with a linear dispersion law, 
when the direction of motion is determined by the value of the spin, the former being able to flip due to magnetic 
impurity scatttering or by coupling with an external field. The resulting intrication between orbital and spin degrees 
of freedom yields a rather rich and complex dynamics with unexpected features, governed by the single parameter 
a = au/v measuring the ratio between the time of fiight and the Larmor period T. 

The most surprising fact is spin-freezing when the width of the initial wave packet becomes quite small: 
narrowing the latter produces stronger and stronger shielding of the spin which can thus keep for ever the memory of 
its initial value. This robustness could be of interest in applications where the spin value carries sensitive information, 
e. g. in spintronics and in quantum spin computation. 

Aside this remarkable fact, the density profile itself shows up a variety of shapes which reflect the competition 
between motion in space and Larmor precession. In the intermediate case where the characteristic time of flight a/v is 
of the order of the precession period, the profile displays a rich structure corresponding to the various quantum paths 
with zero, one, two,. . . precessions: the number of rotations can be directly read by looking at the maxima of the 
density. In one extreme case (a <C 1), one recovers conventional Quantum Random Walk in a space-time continuous 
framework, and its characteristic multimodal (bimodal for J — 1/2) distribution. In the opposite case a ^ 1, the 
particule hardly moves in the mean because of rapid precession, but the width of the packet increases proportionally 
to the time t. 

Another interesting fact, already discussed in the past in the restricted QRW limit, is the role of the phases 
present in the initial spin state: it turns out that the space-density is even only when the initial state is an eigenvalue 
of the transverse spin Jy^ i. e. when the coefficients have the same modulus and a definite phase relationship. 
For an infinitely narrow wave packet, this specific property disappears, and the subsequent motion becomes phase- 
independent: we do not have a simple explanation of this phase symmetry breaking. 
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In all cases, the RMS deviation of the coordinate Ax{t) increases like t, much faster than in classical random 
motion. This universal increase integrates in fact various shapes, from sharp rapidly moving well defined peaks, to a 
standing flattening gaussian distribution. The same can be said about the Shannon entropy, which does not show up 
any cross-over when a varies; it essentially behaves like In Ax, as is often the case, and thus increases oc i at times 
large enough. 

Some of our results are restricted to the J = 1/2 case. Obviously enough, generalization to arbitrary J is 
appealing; in particular, it would be interesting to look at the high- J (quasi-classical limit), especially in order to 
analyse the above mentionned phase symmetry-breaking phenomenon. Furthermore, it would be interesting to check 
whether the above shielding effect is robust against phase decoherence, i. e. to develop simple models incorporating 
coupling to a quantum or classical bath. Work in these directions is in progress. 

Appendix 

We now briefly sketch the methods allowing to obtain the approximate expressions (|50|) and H51|l given in the 
main text. In all cases, the approximations start with the expression of the time-evolution operator U , and preserve 
unitarity. 

Let us begin with the easiest case, namely a ^ 1, i. e. when the spin precesses quickly and when the drift is 
slow. Then, the angle 9 in l|21|l is close to ^. From (|45|l . one readily obtains: 

-1 /--t-oo 

M-, t) - J^^^^H^, J ^ dfce-'=^+^'=^ ,-i.*yiTW^/2 (57) 

and ip-.{x, t) = i'0_|_(— x, t). Taking advantage of the gaussian cut-off, the square-root can be safely expanded when 
— '-^ — <^1; the resulting gaussian integrals are readily evaluated, and one eventually finds the result given in (jSlfl . 

The other case, a ^ 1, (swift drift ans slow precession) is much more involved, due to the fact that the limit 
a ^ is highly singular. From H21|l . the expression of 4'+ is: 

^+(-, t) - (,^)3/4,i/2 y_ dfce-'=^+''^^e-'-^/^^ , (58) 



and one still has \'ip^{x, t)\ = |^+(— a:, t)\. No simple approximation seems possible on such an expression but, since 

2 



one expects that the rapidly moving peak aroud x ~ ^ be slighthy modulated by the precession, it is tempting to 



write an approximation using a convolution integal. For that purpose, we rewrite H58() as follows 



1 /• + 00 

V.+ (x,<) = (3^)3/4^1/2 j dfce"=^e-'= -i'^f^e'^C^-^/^^^) . (59) 



By the convolution theorem, this means that '0+ can be expressed as 

f + 00 



1 p+ca 

^+(-' *) = 23/4,1/4,1/2 j ^ dX'e-ii^-^'-^) F{X', t) , (60) 

where F{X, t) is the Fourier transform: 

F(X, i) ^ — / dfce"=^e'^('=-^'='+"') . (61) 

Note that for a = 0, F{X) reduces to the Dirac function 5{X). Now, we expand the (small) phase factor as 
k — vP+c? ~ k — |fc|(l -t- ^), allowing to write: 



where the function / is defined as: 



vt 
FiX,t) ~ f{X,t)-f{X + —,t) (62) 



f(X, t) ^ — ^k^iX[k^a-vt/(4aXk)] (gg) 



2vr.,o 
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FIG. 7: Kernel F{X, t) defined in ^ for the convolution (see ^). 



This function can be expressed with the Hankcl functions hI ' (2)^^. A somewhat tedious calculation yields 



f{X, t) = -sgn(X) 




(64) 



By (|62|) . this completes the determination of the kernel F{X, t); nevertheless, analysis reveals that the second term 
in (|62l) is quickly negligible, because of the rapid translation of each well-defined peak, so that one nearly always 
has F{X, t) ~ f{X, t). Fig. [7| shows the variation of \F(X, t)\ as a function of the reduced abscissa X = -, at 
short and long times. In all cases, the strong peak in F{X, t) near X — explains the persistence of the two main 
peaks in the density P{x, t), whereas slowly decreasing {^ X^^^^) oscillations are responsible for the secondary small 



r(2) 



peaks located just behind the main one. The minima arise from the zeroes of the Yi Bessel function included in H^ 
which becomes denser and denser as time goes on; this can explain that at large times, the moving peaks for a ^ 1 
eventually grows up (see Fig. 0] upper left). 

I am indebted to C. Caroli, R. Mosseri, P. Ribeiro and J. Vidal for helpful and numerous fruitful discussions. 
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